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Abstract: In this paper we investigate gravitational waves production during 
the Big Bounce phase inspired by the Loop Quantum Cosmology. We consider 
the influence of the holonomy corrections to the equation for tensor modes. We 
C"-- ■ show that they act like additional effective graviton mass, suppressing gravita- 

. . tional waves creation. However, this effects can be treated perturbatively. We 

investigate the simplified model without these corrections and find its exact ana- 
lytical solution. For this model we calculate a spectrum of the gravitational waves 
from the Big Bounce phase. The obtained spectrum decreases to zero for the low 



energy modes. Based on this observation we indicate that this effect can lead 
to the low CMB multipoles suppression and gives a potential way to test Loop 
Quantum Cosmology models. We also consider a scenario with a post-bounce 
inflationary phase. The obtained power spectrum gives qualitative explanation 
of the CMB spectra, including low multipoles suppression. This result is a chal- 
lenge to construct a consistent bounce+inflation model in the Loop Quantum 
Cosmology. 
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1. Introduction 

Gravitational waves seems to be the best tool to explore early universe. In par- 
ticular they can be potentially used to verify quantum cosmological models. This 
idea is based on the fact that gravitational waves produced during the quantum 
epoch can survive frozen on the super-horizontal scales. Then after re-entering 
the horizon they can give imprint on the CMB spectrum, which is observed to- 
day. From the empirical point of view, the most promising is the B spectrum 
of the CMB polarization. This spectrum has its source only in the tensor part 
of perturbation and when observed gives the direct method to investigate relic 
gravitational waves. It is expected that the PLANCK mission could give an op- 
portunity to detect the B polarization 0. However, such predictions base on the 
simple inflationary models. 

One of the most promising approaches to quantize gravity is Loop Quantum 
Gravity . Based on this the theory of the quantum universe, Loop Quantum 
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Cosmology || arose. This theory predicts that the initial singularity state is 
replaced by the quantum Big Bounce ||. In this scenario a universe undergoes 
contraction and then after quantum bounce evolves toward an expanding phase. 
During the the bounce energy density reach maximal finite energy density p c . 
The phenomenological description of the bounce phase can be obtained from the 
modified Friedmann equation 



H< 



-,P 1 



P_ 



;i.i) 



where k = 8nG. 

Investigation of the perturbations in the cosmological models are crucial from 
the point of a large scale structures creation and exploration of the early universe. 
In the Loop Quantum Cosmology this issue has been studied in Ref. [|, || [7| §] . 
In the present paper we consider a particular kind of metric perturbations, the 
gravitational waves. The creation of gravitons in models inspired by the Loop 
Quantum Cosmology were initially studied in Ref. [|, |10| . Then the equation for 



the tensor modes with holonomy and inverse volume corrections has been derived 
in Ref. II]. Recently the equation with holonomy effects has been applied to 



the inflationary phase |L2 



In the present paper we investigate the creation of the gravitational waves 
during the Big Bounce phase. In our considerations we take into account only 
holonomy effects. The inverse volume corrections exhibit the fiducial cell depen- 
dence and are not adequate in the models with the flat background |13|]. The 
production of perturbations in the bouncing cosmologies has been recently stud- 
ied in the different context in Ref. fl4|, [15], [16, 17 1. 



2. Tensor modes with holonomy corrections 

The equation for tensor modes with holonomy corrections has been derived by 



Bojowald and Hossain |ll|. In the source free case these equation takes a form 

(2.1) 

where i = ©, £g> indicate polarization state and 

sin (fl'jk) 



d 2 , /sin2u7£A d , „ 9 , m , 



^- 2 |jf),v 
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(2.2) 
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Here p = a 2 and k = p'/2p. The Hamilton equation for the variable p takes a 
form HIT 



what leads to the 



-/ -^sin2/r)/fc\ 



/sin2/i7fc\ a' 

I 2(i7 / a' 



(2.3) 



(2.4) 



The above equality indicates that the friction term in equation (|2.1| ) holds classical 
form. The effects of holonomies is the additional correction term Tq to the 
classical equation for tensor modes. This factor acts like an additional effective 
graviton mass. 

To define the correction Tq we need to specify a function JX. In general there 
is some freedom of the choice of this function in the power law form. However, it 
has been recently shown that for the flat FRW models the only consistent choice 
is 



1* 



A 
V P 

where A = 2a/3tt7^pi, which is called a /i-scheme. 
Introducing a new variable 



(2.5) 



aha 



ah 6 



u 



and performing the Fourier transform 

d 3 k 



(2.6) 



U T,X 



(27T) 



7 u(t, k)e" 



k-x 



(2.7) 



we can rewrite the equation Q2.1|) to the form 

d 2 



dr 



2 u(t, k) + [k 2 + m 2 eS ]u(r : k) = 



where k 2 = k • k and 



m 



cff 



To . 



(2.8) 



(2.9) 



To calculate this function we must to specify the background dynamics. We 
consider the model with a free scalar field. In this case evolution of the parameter 
p takes a form [T^ 



p = (A + Bt 2 ) 



2\V3 



(2.10) 
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where 



1 3 

A = -/ot?7 2 A and B = -aot?. 
q <f>> 2 * 



Based on definition (12.91) we calculate 



(2.11) 



m 



eff 



k^I (t 2 - f 7 2 A) 
4 (A + Bt 2 f /r 



In the case Tn = we obtain 



(2.12) 



"&(Tg = 0) 



/s: 2 ^ (t 2 - | 7 2 A) 
4 (A + Et 2 ) 573 ' 



(2.13) 



In Fig. D we show the evolution of the effective masses m 2 ff and rn 2 cE {TQ = 0). 
Values of these functions are crucial from the point of view of gravitational waves 
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Figure 1: Evolution of the effective masses m 2 ff and m 2 s (TQ = 0). 

production. Generally, more negative this function is more gravitational waves is 
produced. We see that the presence of correction Tq leads to suppression of the 
gravitons production. However, this effect is relatively weak. 

For the bouncing universe considered we can write a closed system of equa- 
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tions 



du 

dr 

d7r u 




L , - 2 K (t 2 - §7 2 A) 


dr 

dt 
dr 


4 (A + Bt 2 f 3 
+ Bt>f\ 



(2.14) 

u, (2.15) 

(2.16) 



to describe the evolution of the gravitational waves. This system contains the 
holonomy correction to the background dynamics as well as to the correction to 
the equation for tensor modes. Numerical solutions of this system can fully deter- 
minate the classical evolution of the gravitational waves during the Big Bounce 
phase. 

In Fig. ^| we show typical solutions for the function u and hi for the vacuum 
initial condition u ~ \j\fk. It is transparent that the tensor modes are amplified 
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Figure 2: Left : Evolution of the field u. Right : Amplification of the tensor modes 
h during the bounce. 

during the bounce phase. 

In the next section we aim to determinate qualitative and quantitative prop- 
erties of the gravitons produced during the bounce phase. 

3. Toy model 

The goal of this section is to give the introduction to the process of particles 
production during the bounce phase. We show a simple, fully analytical model 
of the gravitational waves creation during the bounce. The model considered 
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does not contain any quantum holonomy corrections to the mode equation. This 
simplification can be however justified. Namely, as it was shown in the previous 
section, the holonomy correction Tq leads only to the perturbative modifications. 
The model considered is based on the two assumptions: 

• we neglect quantum holonomy corrections Tq in the mode equation, 

• we assume matter content in the form 

P=~ 2 - (3-1) 

3.1 Background dynamics 

Plugging condition ([3.1D to the modified Friedmann equation (|1 . 1| ) we find the 
solution 

a(t) = v 7 ! + (t/t ) 2 (3.2) 

where 

'« = 8SV (3 - 3) 

For the further calculations it will be useful to express this solution in terms of 
the conformal time. Performing transformation dr = dt/a we find 

a( T ) = cosh(r/t ). (3.4) 

3.2 Mode functions 



Equation (|2.8|) without the term Tq takes a form 

d 2 , , , r , n a 



, 2 u(r,k) + [k 2 --]u(r,k)=0. (3.5) 

dr z a 



This equation can be obtained from the action 

1 



S t = ^ J dTd 3 x.[u' 2 - 5 ij diud jU - m 2 cS u 2 ] (3.6) 



where 



2 a 



< ff = • (3-7) 

a 



Canonical momenta conjugated to the variable u we obtain from 

6S 



7 =u. (3.1 
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The next step is to quantise this theory. We perform canonical quantisation 
(u, 71") — > (u, 7r) introducing relations of the commutations [u(x, r),7r(y, r)] = 
i5( 3 )(x — y) and [w(x, r), w(y, r)] = [7r(x, r), 7r(y, r)] = 0. The operators w, 7r can 
be decomposed for the Fourier modes 



U{T, X) 



7r(r, x) 



2(2tt) 3 / 2 
1 



2(2vr) 



3/2 



rf 3 k 

d 3 k 



Wk(r)e ik ' x + 4(r)e- tk - x 



7r k (r)e 



ikx | ~i ( \ — ikx 



7r k (r)e~ 



where 



Wk(r) = Ok/ (A;, r) + a f _ k /*(^, r), 
7Tk(r) = a k 0(fc, r) + ol k </*(A;, r), 

and f(k,r)' = g(k,r). 

The equation for the mode function takes a form 

d 2 



dr' 



;f(k J T)+[k 2 +m 2 cS }f(k,T)=0. 



where 



m 



cff 



a 
a 



- - -k 2 



With the use of definition 



i n 



^vrGp c 



and p c 



V3 



16tt 2 7 2 ^ ' 
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and assuming 7 = j M = 0.12738 pOf we obtain 

2.38 



/.'n 



PI 



(3.9) 
(3.10) 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



Now we can find solutions of equation ( |3.11|) . We consider two cases k 2 > /c 2 



and k 2 < k%. 



v_y£LSG Ki J^ Kn 



The solution of equation (|3.11|) takes a form 

f(k, r) = Ae- inr + Be tnr 
where A,BeC and 



n 



U2 _ Ul 

rv thr\ . 



(3.15) 
(3.16) 
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Performing the normalisation, with the use of the Wronskian condition, we 
obtain 



MP - \B\ 2 = 4-. 
11 ' ' 2Q 



We choose advanced modes taking 



.4 



^n 



and 5 = 0, 



what gives 



f(k,r) 



—iCl-r 



g(k,r) = f'{k,r) = -i 



-iClr 



(3.17) 



(3.18) 



(3.19) 
(3.20) 



v^ase K \ Kq 



The solution of equation ( |3.11| ) takes a form 

f(k,r) = Ae~ ClT + Be €lT 
where A,5eC and 



n= xlkl-k 2 . 



(3.21) 



(3.22) 



Performing the normalisation, with the use of the Wronskian condition, we 

obtain 

i 



BA* - AB* 



2Q 



what is fulfilled by 

In this case we obtain 

f(k,r) 



* , „ 1 

A = — == and B 



2VU 



2Vtt' 



2VU 



e Qr + ^-Clr 



g(k,T) = f'(k,r) 
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(3.23) 
(3.24) 

(3.25) 
(3.26) 



3.3 Power spectrum 

The correlation function for the tensor modes takes a form 



(0\K(x,T)h b Mr)\0) = 4 



16ttG 



d 3 k 



a* J (2tt) 3 
(i/c ., . sinfcr 



\f(k,T)\'e 



2„— ik-r 



(3.27) 



where we have defined the power spectrum 



^ ., . 64ttG k 3 . ... .., 

p T (£;,r) = — ^ — /(fc,r)| 2 . 



(3.28) 



Our goal now is to determinate this spectrum on the Hubble (horizon) scales. 
The evolution of the Hubble radius during the bouncing phase considered takes 
a form 



ds 2 = 



i?H — ±o(t) / dr' = ±cl(t)t 
Jo 



(3.29) 



where + denotes expansion and — denotes contraction. We can also express it 
in term of the coordinate time t, obtaining 



.Rh = ±a(t) / 
Jo 



dt! 
'ajt r ) 



±tc 



(t/t ) 2 arcsh(t/t ) 



(3.30) 



This function decreases to zero in the pre-Big Bang phase and increases in the 
post-Big Bang phase. In Fig. || we show evolution of the horizon for the typical 
bouncing cosmologies. We draw also evolution of the arbitrary length scale A oc a. 
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Figure 3: Horizon in the bouncing cosmologies. 
We see that during the bouncing phase all length scales finally cross the horizon. 
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Quantum fluctuations of the length A < Rn behave like a virtual particles, when 
fluctuations of the length A > i?n become classical excitations. 
The condition A = Rn indicate 

2na 2na 2tt 



k 



(3.31) 



Plugging it to definition ( |3.28| ) and using the mode functions ( |3.19| ), ( |3.25| ) we 
determinate the tensor power spectrum on the horizon scales 



2 / \ 3 
16 / fc \ I k 
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(3.32) 
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We plot this spectrum in Fig. £|. We see that for k < kn the spectrum is charac- 
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Figure 4: Spectrum of the tensor perturbations on the horizon. 

terised by the bump. For k — > the spectrum is decreasing to zero. However, for 
k — > kn the spectrum has a pole Vt — > oo. This pole is however not physical. To 
show it we consider the correlation function 



where 



(0\hl(x,T)h b a (x,T)\0) 



\f(k,r)\ 2 = Q(k-k ) 



32G 



it a* 



dkk 2 \f(k,r) 



+ B(k - k) 
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ch 



2x1 k\ - k 2 T 



(3.33) 
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The correlation function is a physical quantity and can indicate which features 
of the spectrum Vt are real divergences. As we can see, near the pole k = k* the 
correlation function is determined by the integrals 



rv rv I rbi 



.2 



Ix =-- / dk = \ x = T~ ( = ~^~ x ^ x2 ~ 1 — In x + a/x 2 — 1 

a/A; 2 - fcjj I «o 



(3.34) 



and 



dk — . = < x 



fcl _kl 



-xy/l — x 2 — arcsinx 



(3.35) 



Both of them are finite for the x = 1 and give a finite contribution to the definition 

of the correlation function. 

For k — >• oo the power spectrum exhibit a UV divergence. In fact, the region 

k > k traces the Planck scales and a different approach should be used to 

determinate properties of the quantum fluctuations. Another important issue is 

to indicate which of the modes survive frozen above the Hubble radius. In fact, 

not all of them but only these with k < k$. It can be seen from the asymptotic 

solutions of equation ( |3.5|) . Namely, in the regime k 2 <^C \a /a\ = kg we have 

solution 

f T dx 

Here we have a constant contribution A^ which freezes the amplitude of the 
gravitational waves. In this regime the spectrum does not change during the 
evolution above the horizon and will be the same when re-enters the horizon. In 
the regime k 2 ^> \a /a\ = k\ we have decaying solutions 



p ±ikr 



(3.37) 



These modes does not lead to the classical fluctuations when re-enter horizon. 

Summing up, the spectrum Vt on the second branch of the horizon is free 
form the UV divergent part for the k > ho- The only contribution to this spectrum 
comes from the bump for the k < ko- 

3.4 Parameter f2 gw 

To describe the spectrum of gravitons it is common to use the parameter 



p* dv 
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where p gw is the energy density of gravitational waves and p* is present critical 
energy density Our goal in this section is to calculate the function f2 gw (z/) for 
the gravitons produced during the Big Bounce phase. 

We consider the creation of the gravitons during the transition from some 
initial to final states. The initial vacuum state |0i n ) is determined by ak|0; n ) = 0, 
where a^ is the initial annihilation operator for Tj. The relation between an- 
nihilation and creation operators for the initial and final states is given by the 
Bogoliubov transformation 



k J 



6 k = B + (k)a k + B„(k)*al 
bt = B + {k)*al + S_(/t)a_ k 



(3.39) 
(3.40) 



where \B + \ 2 — |-B_| 2 = 1. Because we are working in the Heisenberg description 
the vacuum state does not change during the evolution. It results that &k|0m) = 
.B_(/c)*a_ k |0i n ) is different from zero when B-(k)* is a nonzero function. This 
means that in the final state graviton field considered is no more in the vacuum 
state without particles. The number of produced particles in the final state is 
given by 

% = -(0in| fc k fck + &liA-i 

The energy density of gravitons is given by 

Airuj 2 duj 



\Q in ) = \B4k)\ 2 . 



(3.41) 



dp. 



S\v 



2-fku- 



(2tTc) E 



\B-(k)\ 2 . 



(3.42) 



where we used definition ( |3.41| ). The expression for the parameter f2 gw defined 
by ( |3.38| ) takes now the form 



^gw(^) 



where 



and 



Vtn ■ u^ ■ sinh 2 



ftn 



he 16ir 2 



c 4 p* 



(k/ko) 
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(3.43) 



(3.44) 



(3.45) 



For the model with the inflationary phase we have — ~ 



10 . Another value which 



must be specified is the duration of the bounce. We assume that n 
and the r f = 20 Zpi- 
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Figure 5: Function J7 gw (^) for the r; = —20 Zpi and the Tf = 20 lp\. 

In Fig. |5| we show the function Q gw (^) for the setup considered. The spectrum 
is characterised by the bump with the maximal frequency 



7-10- 14 [Hz]. 



(3.46) 



The presented result approves our previous consideration. Namely, there 
is no production of the gravitons with frequencies above some ^ max - The only 
contribution to the spectrum comes from the bump. 



4. Tensor modes from the pre-bounce phase 

The bouncing solution ( |2.1U| ) in the limit t — > ±oo gives 

a(t) oc It] 1 / 3 oc M 1 / 2 . 



(4.1) 



In this regime the correction Tq vanishes and the expression for the effective mass 
simplifies to 

2 a 11 

meff= ~T = 4^- 
The equation for the mode functions takes a form 



(4.2) 



dr 2 



f(k,r) + 



k 2 



1 1 

4^2 



/(*,r) = 0, 



The normalised solution of this equation has a form 

A' ,^. T (i), 



f(k,r) 



V^kH ( L> (-rk) 



(4.3) 



(4.4) 
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where 



M 



,J7r/4 



(4.5) 



With the use of definitions ( |3.28| ) and 
pre-bounce phase in the form 



we obtain the power spectrum for the 



V T {k) 



12, 




r(l) 



7T 



H^(2k)\' 



k 

k# 



for k —>■ 0, 



(4.6) 



where k# is some constant. 

It is important to note that this part of the spectrum does not depend on 
the quantum gravitational effects. It is in opposite to the predictions from the 
inflationary models, where low energy modes come from a high energy region. 
Here, low energy modes are produced in the low energy pre-Big Bang state. So, 
the present largest scale structures have their origin in the semi-classical pre-Big 
Bang phase rather than in the deep quantum regime. 

5. Suppressing low CMB multipoles with a bounce+ inflation 
scenario 

In this section we consider scenario with an inflationary phase taking place after 
the bounce. In Fig. ^| we show the resulting evolution of the horizon. We see that 




Bounce 



Inflation 



Figure 6: Evolution of the horizon in the bounce+infiation model. 

the horizon is firstly crossed by the modes of the length A > A*. These modes 
lives frozen on the super-horizontal scales and re-enter the horizon in the post-Big 
Bang phase. Additionally they are not sensitive for the quantum gravitational 
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effects close to the bounce. As it was shown in the previous section they leads to 
the power spectrum in the form 



V T (k) oc k 3 for k < k* 



(5.1) 



These fluctuations can lead directly to the CMB fluctuations on the large angular 
scales. Modes of the length A < A* also cross the horizon in the pre-Big Bang 
phase but re-enter shortly after bounce. They decay shortly after this and be- 
come a quantum fluctuations. Then, during the inflationary phase they cross the 
horizon again. The spectrum of these fluctuations will be nearly flat as predicted 
by the inflationary models 



V T (k) = A for /c> £;*. 



(5.2) 



The loop holonomy corrections to the inflationary spectrum were studied in Ref. 

Finally, we should obtain the spectrum of the tensor perturbations in the 
form shown in Fig. |7|. 



V = const 
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Figure 7: The schematic picture of the power spectrum expected for the 
bounce+inflation model. 

To the describe power spectrum in a continuous way we propose the phe- 
nomenological formula 

k 



V(k) =A[- 



!+(£ 



3' 



which interpolates the two regimes considered. We plot this function in Fig. R 
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Figure 8: Phenomenological spectrum for the bounce+inflation model with A = 1 
and k* = 1. 



The spectrum in the same form should be also obtained for the scalar modes 
of fluctuation. In this context similar models were studied in Ref . [f2^, |24| . They 
have shown that a transitional regime in the spectrum can have form of the 
oscillations, as it was shown in Fig. [?]. However, in the cited papers evolution 
was obtained form the specific dynamics of the single inflationary field. Namely, 
no quantum gravitational corrections was used to obtain the bounce phase. 

An interesting property of the obtained spectra is damping of the low mul- 
tipoles in the CMB spectra. Such an effect is in fact observed p5| . This gives 
potentially the possibility to test bouncing cosmologies. If the bounce is present 
we should observe low multipoles suppression. However, this effect can have also 
other origins and on the present stage of observations it is not possible to indicate 
that it is truly the remnant of the Big Bounce. 

The important feature of the presented scenario is that none of the primordial 
perturbations produced in the deep quantum gravitational regime gives seeds to 
the structures formation. All perturbations come from either the pre-Big Bang 
semi-classical phase or from the post-Big Bang inflationary phase. It is therefore 
hard to investigate observationally the deep quantum regime. However, we can 
potentially observe the classical pre-Big Bounce branch which is the result of the 
quantum gravitational effects. 



6. Summary 

In this paper we have considered the gravitational waves creation during the 
Big Bounce phase inspired by the Loop Quantum Cosmology. We have studied 
effects of the holonomy corrections to the equation for tensor modes. We have 
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shown that they lead to the suppression of the gravitons production. However, 
this effects is not dominant and can be treated perturbatively. In fact, to obtain 
qualitative results it is justified to neglect this contribution. Based on the above 
studies we have solved the simplified model of the gravitons production during the 
bounce phase. We have derived values the power spectrum Vt and the parameter 
Q gw . The obtained spectrum for the bounce has a form of the bump. It decreases 
to zero for the energies tending to zero and for some high energy scale. This is 
a typical property of the bouncing cosmologies. The similar behaviour is also 
expected for the scalar perturbations. The perturbations form the pre-Big Bang 
phase gives a direct imprint for the large scale structures. In particular they can 
lead to the suppression of the low CMB multipoles. However, the bounce phase 
itself does not lead to observed features of the CMB spectra. The post-bounce 
inflationary phase is required. Such a model with the bounce and the following 
inflationary phase is not yet constructed in the Loop Quantum Cosmology. It is 
therefore a challenge to construct it in a consistent way. 
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